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ABSTRACT 


A method for estimating extinction coefficients in the near-infrared and infrared 
wavelengths from satellite measurements in the visible was developed. Five tests were 
devised to examine the limits and sensitivity of the model. The first test studied the 
error inherent in the retrieval of the parameters which are needed to describe the 
distribution of atmospheric particles and are direct inputs for the calculation of 
extinction at 3.75 and 11.00 tum. Also studied were errors associated with uncertainties 
in the extinction values, uncertainties in relative hunudity values, deviations of particle 
size distribution from the model and effects of high winds on the aerosol distribution. 
Results indicate that the biggest error results when wind generated aerosols change the 
particle size distribution especially at radii larger than 2 um. The error reaches 329^ 
for prediction at 11.00 um at 60% relative humidity. The smallest error, less than 7% 
for all variations, is associated with the retrieval technique itself. Errors up to 25% in 
the measured satellite extinction coefficients lead to errors 01 ظا‎ to 2320 1 E 
estimated values for both the marine and rural models. Results indicate that the rural 
model at high values of relative humidity is affected the most with an error of 31% ata 
RH of 95% at 11.00 tm. Negative deviations in the marine particle size distribution 
give rise to large errors for 入 = 0 um. For 20% deviation, the error Cal Ts 


high as 41%. The error decreases accordingly as the percent deviation 1s reduced. 
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I. INTRODUCTION 


A. MOTIVATION AND BACKGROUND 

Electro-optical surveillance, guidance and weapon systems are affected by 
oceanographic and meteorological parameters (Mitchell er al., 1982). Environmental 
features such as clouds, rain, relative humidity, haze, dust, aerosols and precipitation 
have been shown to have significant effects on the performance of these svstems. 
Mitchell's study points out that systems such as the REWSON optical missile warning 
set used by helicopters, medium range air-to-surface missiles, the AEGIS combat 
weapon system (a fullv automatic surface-to-air system) and high energy lasers have 
been found to be susceptible to the presence of aerosols in the atmosphere. These 
environmental parameters are not always considered in evaluating the performance of 
these weapon systems. This is because either the environmental information is not 
available, or the scale of the phenomena is so small that little confidence can be placed 
in the accuracy of its specification. In order to improve the evaluation of the 
performance of these systems the applicable parameters should be included for each 
weapon system. 

Models such as the LOWTRAN 6 used for the evaluation of performance of 
certain missile systems can handle properly the absorption and scattering of radiation 
by water vapor for the prediction of transmission losses but improvement could be 
done in the initialization to account for a better representation of aerosols (Kneizys et 
al., 1983). A joint effort bv the Naval Research Laboratory, the Naval Postgraduate 
School and the Naval Ocean Svstems Center is presently underway to develop a Navv 
Ocean Vertical Aerosol model C(NOVAM) for inclusion into a future version of 
LOWTRAN. A different approach involving small modifications of LOWTRAN 6 
could be done that would allow the insertion of an extinction coefficient value at the 
infrared wavelengths rather than the actual parameters needed to calculate the particle 
size distribution. This extinction coefficient value would be obtained from sateilite 
estimates of optical depth and would take into account the aerosols present at the 
target location rather than at the point from where the weapon system is launched. 
This 1s important because some weapon systems rely on target signature in determining 


system performance. Target signature is related to the target-to-background 


temperature difference, which, when propagated through the atmosphere, is modified 
by the transmittance and thus by the presence of aerosols. 

Since there are weapon systems that operate at infrared, near-infrared and visible 
wavelengths, evaluation of the effect of environmental parameters must be performed 
at more than one wavelength. Estimating the extinction of the atmosphere (related to 
parameters such as visibility) must be done at the visible and near-infrared wavelengths 
and extrapolated to infrared wavelengths. This infrared extinction may be obtained 
from satellite measurements, but at the infrared wavelengths the power received from 
scattered radiation is small compared to visible wavelengths. Also, due to the presence 
of ambiguities from a variable background and water vapor, it is more practical to 
obtain measurements at visible wavelengths and estimate the extinction at the infrared 
wavelengths. 

Studies by Gerber (1985) show that it is possible to use visible light scattered by 
aerosols to give useful estimates of extinction at infrared wavelengths. He showed that 
measurements using a combination of a two channel integrating nephelometer, which 
operates at one wavelength over two angular ranges, and an axial-scatter sensor gave 
the best accuracy (~ + 109^) for scaling over the 1.06-10.06 ttm range for extinction by 
marine aerosols. However. his studies require direct measurement of the scattered light 
by the above instruments. The object of this thesis is to obtain an accurate 
representation of the size distribution of the aerosols in the marine atmospheric 
boundary laver from satellite measurements of radiance at visible wavelengths. The 
size distribution will then be used to obtain estimates of extinction coefficients at near- 


infrared and infrared wavelengths. 


B. ORGANIZATION 

The thesis is composed of five chapters. These include an introduction, theory, 
procedure, results and discussion, and a summary of conclusions and 
recommendations. Chapter II will discuss the theory of radiative scattering by 
atmospheric aerosols and the basic assumptions used in computing the scattering and 
extinction coefficients. A brief description of the atmospheric aerosol models is also 
included. Chapter III describes the model for estimating extinction coefficients. at 
near-infrared and infrared wavelengths. It discusses the assumptions dealing with the 
vertical distribution of atmospheric particles in the marine boundary laver. Also 


included is the retrieval of the parameters required to describe the particle msnm 


10 


distribution of the aerosois and the caiculation of extinction at 3.75 and 11.00 ۰ 
Chapter III also discusses the validity of the numericai techniques used. Chapter IV 
discusses the procedure and results of five tests developed to evaluate the sensitivity 


and limitations of the models. A summary of the results and recommendations are 
stated in Chapter V. 
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IIl. THEORY 


A. EXTINCTION AND SCATTERING COEFFICIENTS 

Radiation can be scattered by air molecules, aerosol particles which are 
suspended in the air and by the earth’s surface. Scattering by particles which are much 
smaller than the incident wavelength, such as molecules, is described by Ravleigh 
theorv. Therefore, in the presence of gas molecules (~ 1074 jm) and small particles 
(70.05 um), radiation undergoes Rayleigh scattering at visible wavelengths. For 
particles whose size are comparable to or larger than the wavelength, the scattering is 
referred to as Mie scattering. Thus, in the presence of cloud droplets (^ 10 nm) 
radiation undergoes Mie scattering at the infrared wavelengths, while in the presence of 
aerosols ( ^ 0.1-1.0 utm) radiation undergoes Mie scattering at the visible wavelengths. 
In a cloud-free environment away from large temperature or pressure gradients the 
molecular constituents do not vary in the horizontal appreciably and thus Ravleigh 
scattering does not vary significantly in the horizontal. However, Mie scattering 
depends on the size and number of aerosol particles. These partcles are primarily 
generated bv sea sprav in a marine environment and are composed primarilv of dust- 
like particles over the continents. Advection of these rural aerosols would result in 
large horizontal variations of the particle size distribution. 

Scattering of light by aerosols is characterized by the extinction coefficient, 


ext 
pss 1s given bv the sum of the extinction by scattering and absorption. 


pe E Dua T B bs (220 


Absorption of radiation by aerosols 1s described by the complex part of the index of 
refraction. For marine particles the complex part of the index of refraction is less than 
107* at wavelengths less than 1 um, and therefore absorption bv marine aerosols is very 
small (Shettle and Fenn, 1979). For visible and near infrared wavelengths, Pest for 


marine aerosols is approximately equal to B For higher wavelengths and for rural 


scat’ 
cases, the complex part of the index of refraction 1s of order 1072 so both terms must 
be included in the calculation of extinction at 3.75 and 11.00 pnm. Gases in the 


atmosphere, such as CO, and water vapor, absorb radiation. However, at the 


wavelengths of interest (0.63, 0.86. 3.75 and 11.00 pum) the absorption by these 
constituents is small and can be neglected in our calculations. 
Given the size and distribution of these particies, the extinction coefficient (in 


units of per length) is defined as: 

Bost * of TO. MDENEOIdD, (2.2) 
where Q,,, is the extinction efficiency of a particle with radius r and complex index of 
refraction m, and dN/dr describes the size distribution of the particles. As seen in Eqn. 
2.2, the extinction coefficient depends on three terms: the cross-sectional area (mr~), the 
extinction efficiency (Q.,,) and the size distribution (dN/dr). The relationship of r, 
ېے ډ)‎ and dN/dr for a typical visual case is presented in Fig. 2.1. Panel b of Fig. 2.1 
shows a particle size distribution for a relative humidity (RH) value of 90° composed 
of two modes. The modes correspond to small (mode 1) and large (mode 2) radii 
particles. The nature of each mode is described in detail in l. B. This example shows 
that when r is small, rr” and Qey, are small, while dN/dr is large: when r is large, nr 
is large and dN/dr and Q,,, 


bounded region of particle sizes. Low values of Q 


become small. The result is that f) 


eum bound DR - 


at the large radii. For visible wavelengths and a 


ext 18 affected by a 
at the small radu 
while a decrease in dN/dr bounds f... 
marine aerosol size distribution model at 80% relative humidity,. studies done by 
Durkee (1987) include a particle range from approximately 0.5 to 5.0 um. To account 
for variations in wavelength and relative humidity this studv incorporates a radi range 
from 0.01 to 50.0 um in the extinction coefficient calculation. 

The extinction efficiency function is dependent on wavelength and relative 
humidity since composition changes with increased condensed water. As wavelength 


increases Q.,, shifts to the right, weighing larger radius particles more than smaller 


Ext 
particles as depicted by Figs. 2.1 and 2.2. Since there are fewer particles at the larger 
radii the net effect of increasing wavelength is a decrease in the extinction coefficient, 
Pext Fig. 2.3 shows the effect of relative humidity on the size distribution of marine 
particles. Notice the increase in the number of large particles for higher relative 
humidity values on Fig. 2.3. Since the Q.,, function at 3.75 am weighs the larger 
particles (greater than the mode radius) more than smaller particles, and there are more 
particles at the larger radii, the net effect is an increase in the cumulative extinction for 


higher relative humidity values. 
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Fig. 2.1 a. Extinction efficiency for } = 0.63 um and RH = 90% for mode 1 (solid) 
and mode 2 (dashed). b. Model distribution of marine particles at RH = 90%. Sum 


SO & 2 (solid), mode | (dashed), mode 2 (dot). c. Cumulative extinction from 
enm 2. 
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Fig. 2.2 Same as Figure 2.1 but for 3.75 }tm wavelength. Note change of scale for c. 


۱3 


epar eco Hi”) 





jp 102 10 


RADIS ir. 


Fig. 2.3. Variation of marine model size distribution of Shettle and Fenn (1979) with 
relative humidity. 
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D. ATMOSPHERIC AEROSOL MODELS سرت‎ EER 
1. Introduction 

Atmospheric aerosol particles in the atmosphere vary greatly in their 
concentration, size and composition, and consequently in their effects on optical and 
infrared radiation. Several tropospheric aerosol models have been developed by Shettle 
and Fenn (1979) which include the dependence of atmospheric aerosols on relative 
humidity. 

The size distribution for the different aerosol models used here are represented 


bv the sum of two log-normal distributions: 


2 N. od ao. 
iz1Mn(10) r e 7 27 


where N is the cumulative number density of particles of radius r; © is the standard 
deviation; r, N, are the mode radius and the total number density. This form of the 
distribution. function represents the multimodal nature of the atmospheric. aerosols 
Which will be discussed next. This study will begin with a bimodal size distribution 
since it has been shown to be generally adequate to characterize the gross features of 
most aerosol distributions (Whitby and Cantrell, 1975). A third mode will be added to 
allow for wind generated, larger than 2 {tm mode radi, particles. The rural model is 
comprised of a small rural aerosol mode (mode 1) and a large rural aerosol mode 
(mode 2). Ihe marine model consists of a small rural aerosol mode (mode 1) and an 
oceanic mode (mode 2). 
2. Rural and Marine Aerosol Models 

The rural model is intended to represent the aerosol under conditions where it 
is not directly influenced by urban and/or industrial aerosol sources. The aerosols are 
assumed to be composed of a mixture of 70 percent of water soluble substance and 30 
percent dust-like aerosols. The marine aerosol model is composed of a sea-salt mode 
and a continental mode which was assumed to be identical to the rural aerosol with the 
excepuon that the very large particles were eliminated since they will be eventually lost 
due to fallout as the air masses move across the oceans. The relative proportions of: 


aerosols of oceanic or continental origins will vary, particularly in coastal regions. To 
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account for these variations. the model permits adjustment of the relative contribution 
by the oceanic and continental modes. The contribution from the marine size 
distribution used in this study varied from Zero to 6% of the total number of particles 
of oceanic origin while the large rural aerosols varied from zero to 0.025%. The small 
radii particles, mode 1, make up the rest and were changed accordingly. 

As relative humidity increases, water vapor condenses out of the atmosphere 
onto the particulates suspended in the atmosphere. This condensed water increases the 
size of the aerosols and changes their composition and their effective refractive index. 
The resulting effect of the aerosols on the absorption and scattering of light will be 
modified correspondingly. Table | shows how the size of a particle changes with 
relative humidity for both the rural and oceanic models. As relative humidity 
increases, the mode radii increases for all modes and thus the number of particles at a 
given radius increase with the larger radii increasing the most. as seen previously in 
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TABLE 1 


CHANGES IN AVERAGE MODE RADIUS (rl AND r2) AS A 
FUNCTION OF RELATIVE KUIDEN 


RELATIVE RURAL MARINE 
HUMIDITY r1( um) 12 ) 11۳۱ ( r2( um) 


09 22:10 . 4300 0. 0 
50% . 02748 NOE 
70% . 02846 . 4571 . 2041 
30% . 03274 .5477 . 3180 
90% . 03884 . 6462 . 3803 
95% . 04238 . 7078 . 4606 
98% .04751 . 9728 . 6024 
99% . 05215 . 7555 . 7505 
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The marine model was modified to allow the addition of a wind: driven third 


mode to the size distribution. As shown in Fig. 2.4 this third class of aerosol consists 


18 


of the largest nuciei originating from the sea surface. It has a mode radius of 2 um 
and the amplitude varies with wind speed (Gathman, 1983). These large particles can 
be important in the propagation of infrared radiation near the sea surface under a 


condition of high winds. 
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Ill. PROCEDURE 


A. LOGIC 

The main goal of this study is to be able to estimate extinction coefficients at 
near-infrared and infrared wavelengths from satellite information at shorter 
wavelengths. Channel 1 (0.55 to .68 um) and channel 2 (0.70 to 1.00 um) of the 
NOAA-7 Advanced Very High Resolution. Radiometer (AVHRR) return image 
information from which a value of radiance is obtained. Since radiance can be shown 
to be a nearly linear function of aerosol optical depth (Durkee er al., 1986), the optical 
depth due to the presence of aerosols can be calculated as discussed in IIT. B. 


The optical depth is calculated from integration of the extinction coefficient: 


t= |B.,,dz (3.1) 


Above the boundary laver. as shown in Fig. 3.1 (Fairall and Davidson, 1986), the 
volume of aerosols decreases due to a decline in relative humidity, and therefore the 
extinction coefficient, will decrease as discussed in Chapter II. A. Thus, the extinction 
above the top of the marine atmospheric boundary layer (MABL) can be considered 
negligible since the number of particles above this height is significantly reduced. The 


integration is thus simplified and D... can be approximated from: 


ext 


Dee (3.2) 


where Z 1s the height of the top of the MABL. Therefore, 1 can be estimated from 
satelhte estimates of t and estimates of Z obtained from models such as the Fleet 
Numerical Oceanography Center's Navy Operational Local Analvsis and Prediction 
System (NOLAPS) (Burk and Thompson, 1982). Later in this chapter it will be shown 
that this extinction value is the input required in the model in order to retrieve a 
particle size distribution which is used to calculate the extinction at other wavelengths. 

A simple technique for estimating extinction coefficient at longer wavelengths 
was attempted by plotting extinction coefficient values in the visible and near-infrared. 


Fig. 3.2 shows the relationship between extinction coefficients at 0.63 and 3.75 um for 
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three values of the contribution of the second mode, N2. Each distribution case is 
varied by changing reiative humidity. The total number of particles, TN, was set at 
4,000 cm3, relative humidity, RH, varied from 0 to 99%, and the contribution of the 
second mode, N2, was set at 0.0, 0.000125 and 0.000250 for each line respectively. In 
order to find extinction values at 3.75 and 11.00 um an infinite number of plots similar 
to Fig. 3.2 would have to be generated, one for each distribution shape. A simple 
prediction using one wavelength is therefore impossible. This result 1s consistent with 
Gerber (1985) who showed for various atmospheric simulations that no correlation 
exists between extinction at visible and near-infra red wavelengths when scaling with 
only one wavelength. 

We have shown that the extinction value depends on the number of particles 
present and on their distribution. Thus, a method using satellite-detected radiance at 
two wavelengths is proposed for estimating extinction at the 3.75 and 11.00 um 
wavelengths. Approximations made in order to convert radiance measurements from 
satellite data to an extinction value are discussed in [II. B. From an input of 
extinction values at 0.65 and 0.86 uim the model can retrieve the total number of 
particles (TN) and the percentage of the second mode (N2). With knowledge of these 
parameters, which are needed to describe the particle size distribution. and an effective 
humidity value for the boundary laver, the model can calculate extinction at 3.75 and 


11.00 uim as discussed in III. C. 


B. CALCULATION OF OPTICAL DEPTH AND EXTINCTION COEFFICIENT 

Channel 1 (0.63 um) and channel 2 (0.86 um) of the NOAA-7 AVHRR satellite 
returns image information in the form of digital brightness counts. A digital number 
(which is proportional to the amount of radiation received at a satellite) is assigned to 
each image element. This inserts some uncertainty since only integers are used in the 
digitizing process. Fig. 3.3 shows the conversion between digital counts and physical 
units for the AVHRR (Durkee, 1984). 

The radiance at satellite altitude (Lc) consists of contributions from the following 


terms: 
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一 一 Na 0 
---- N2 0.000125 
N2 0.000230 


Cumulative Extinction (km ') at 0.6328,ım 
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Fig. 3.2 Estimation of extinction coefficient. at 3.75 pm from 0.63 pm for three 
different particle size distribution, described by N2. 
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Ea Conversion of digital counts to radiance for the AVHRR sensors (Durkee, 


Where Ly; is the ‘water leaving’ radiance caused by subsurface reflectance, Lg is the 
contribution from the specular reflection off the sea surface, Lp is the path-added 
radiance due to Rayleigh scattering in the atmosphere, La is the path-added radiance 
due to Mie scattering in the marine boundary layer and T is the atmospheric 
transmittance. 

Ron wavelengths in the red-visible to near-mirared, Ly generally will be a small 
contribution to upwelling radiance. Lg, also referred to as the radiance added by 
'sunglint, can be estimated from knowledge of the sun-earth-satellite geometry and 
surface roughness. Lg is negligible under all but certain geometries. Lp is obtained 
from knowledge of the Rayleigh optical depth (which can be estimated from profiles of 
pressure and temperature) and scene geometry (McCartnev, 1976). 

In the single scattering approximation we consider the radiation which is 
scattered only once by marine aerosols. If we assume that the upward intensity at the 
bottom of the atmosphere is zero and for atmospheres with small optical depth, it. can 
be shown that radiance is a linear function of optical depth (Durkee ez al., 1986).. The 
optical depth due to the presence of aerosols can then be estimated from the aerosol 
radiance. Eqn. 3.2 can then be used to calculate the extinction coefficient required for 
input into the model. For the purpose of this study it is assumed that the extinction 
values at 0.63 and 0.86 ım are available from estimates of’ optical depth (t) and 


boundarv laver depth (Z). 


C. ESTIMATION OF [IR EXTINCTION COEFFICIENT 

Fig. 3.4 shows a basic outline of the model for estimating extinction at near- 
infrared and infrared wavelengths. This model can be divided into two sections. In the 
first part of the model. input of the aerosol model type (marine or rural) and the 
effective relative humidity for the MABL defines the location of the average mode 
radius, rl and r2. Input of the extinction coefficient at 0.63 and 0.86 Jim will determine 
the magnitude of the retrieved total number of particles, IN, and percentage 
distribution of the second mode, N2. These are the parameters necessary to describe 
the distribution of the aerosols particles. The second part of the model uses these 
values of TN and N2 and an effective relative humidity value for the boundary layer to 


estimate the extinction coefficient at 3.75 and 11.00 jtm. 
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11.00 um 


Fig. 3.4 Outline for the estimation of extinction coefficients at NIR and IR 
wavelengths. 
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l. Retrievai of TN ana Nz 

This portion of the mode! is used to calculate the total number of particles 
and the contribution of the second mode needed to describe the aerosol distribution 
using a bilinear interpolation. The model's application is limited by the total number 
of particles and the variations allowed. For this study, in the marine case, TN was 
allowed to vary from 1,000 to 7,000 cm"? and the second mode contribution, N2, 
varied from 0.0 to 0.06. For the rural case. TN was allowed to vary from 10,000 to 
20,000 cm”? and the second mode contribution, N2, varied from 0.0 to 0.000250. 
These are the current limitations of the model but the model can be easily expanded to 
allow for greater variations. The data files are comprised of 7 X 7 pem matrices for 
marine and rural cases at 0.63 and 0.86 um, and for RH = 0, 50, 70, 80, 90, 95, 98 and 
99%. These matrices were obtained bv running the second portion of the model for 
the following values of TN and N2. For the marine case, N2 values used were 0.0, 
0.01, 0.02, 0.03, 0.04, 0.05 and 0.06. For the rural case N2 values used were 0.0, 
0.000050, 0.000100, 0.000125. 0.000150, 0.000200, 0.000250. In order to run the model, 
an effective humiditv for the boundarv laver and values of extinction at 0.63 and 0.86 
jim must be entered. A subroutine in the model, using a Hermite interpolation, allows 


calculation of a 7 X 7 matrix for any integer value of relative humidity less than 


ext 
100% not m the data files. 

As discussed in IH. A., Doer is a function of TN, N2, RH and X. An infinite 
number of contributions of TN and N2 will give rise to a predetermined extinction 
coefficient at a given wavelength. But onlv one TN and N2 conibination will give rise 
to a given B 


at A, and a given Pey, at A, Therefore, we must generate lines 


eXt ext 
describing these combinations at each wavelength and solve for the intersection of 
these two lines which will give the value of TN and N2 associated with the input Des. - 
at the two visible wavelengths. 

The following description shows how the model mathematically solves for the 
two parameters, TN and N2, needed to describe the aerosol distribution in question. 
For a given extinction value, a corresponding TN value is calculated at each N2. This 
generates a set of seven values of TN, one at each N2 value, for each A. The 
intersection between the curves formed by these sets (Figs. 3.5 and 3.6) represents the 
best approximation of TN and N2 which describes the particle distribution 
corresponding to the extinction coefficient value observed at 0.63 and 0.86 um.. This 


single value of TN and N2 is obtained by a bilinear interpolation. The bilinear 


28 


interpolation consists of obtaining an equation for each line using a point-form for the 
nearest point on each side of the intersection. By soiving tne two equations the 
intersection of the two curves is determined and therefore an estimate of TN and N2 is 
obtained. 

2. Calculation of li at 3.75 and 11.00 pm 


This portion of the model calculates an extinction coefficient value, B at 


ext? 
near-infrared and infrared wavelengths, from knowledge of the particle size distribution 
and an effective value of relative humiditv for the marine boundary layer. Data files 
containing mode radii for marine and rural cases are supplied. These values have been 
interpolated using a Hermite interpolation for integer values of relative humidity 
ranging from 0 to 99%. Also included are data files containing values of index of 
refraction for marine and rural cases. These values have also been interpolated using a 
Hermite interpolation for integer values of relative humiditv ranging from 0 to 99%. 
Uninterpolated values for mode radu and index of refracton at RH = 0, 50, 70, 80. 90, 
eo. 97, 98 and 99% were obtained from Shettle and Fenn (1979). The value of TN 
and N2 obtained from the previous portion of the model and an effective relative 
humidity, RH, for the marine boundary layer are entered into the model. Using Eqns. 


2.3 and 2.2 the model calculates D t 3.75 and 11.00 um. 


ext 2 


-D. VALIDITY OF NUMERICAL TECHNIQUES 
Values for the average radius of each mode, ri and r2, and the real and 
imaginary parts of the index of refraction for the marine and rural cases are provided 
for only certain relative humidity values as discussed in the previous section. Figs. 3.7 
and 3.8 are a plot of the average mode radius as a function of relauve hunidity and 
show a monotonic behavior for both ri and r2. Due to the (lat section followed bv a 
rapid nonlinear portion for relative humidity values greater than 70%, a Hermite cubic 
numerical interpolation was used to generate values of rl and r2 for relative humidity 
values ranging from 0 to 99%. Figures 3.9 and 3.10 show a similar monotonic trend 
for the scattering extinction for rural and marine cases as a function of relative 
humidity. Based on this result, matrices needed for extinction values during the 
retrieval portion of the model are generated by a Hermitic interpolation. The real and 
imaginary portions of the index of refraction also follow this behavior and a Hermitic 
interpolation was therefore employed. The Hermitic cubic function in this program 


was developed bv the Lawrence Livermore National Laboratory (Fritsch, 1982). Such 
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an interpolating screme 1s moere reasonable than a cubic spline if the data contain both 
“steep” and “flat” secticns as shown in Figs. 5.7 - 3.10. 

As discussed in the previous section, a bilinear interpolation was used to solve for 
the intersection of the two lines and retrieve the parameters, TN and N2, needed to 
describe the particle size distribution of aerosols. This technique was appropriate for 
this case, as seen in Figs. 3.5 and 3.6, since the equation of the lines are generated with 
points in the immediate vicinity of the intersection and a linear approximation between 
such short distances can be made. Fig. 3.5 shows that if the points were separated a 
greater distance a linear approximation would not hold true for the marine case since 
over the interval used for N2 (0.0 to 0.06) the curve is not linear. However, it would 
still hold true for the rural case since the lines are almost linear as seen in Fig. 3.6. For 
the rural case the values immediately following the intersection (Al-A4) were used to 
find the intersection point. A set of values (B1-B4) further away from the intersection 
Were used and compared with the first set. Prediction errors for extinction at 3.75 and 
11.00 pm were calculated. The results shown in Table 2 show excellent agreement 
since the errors are less than 1% for both infrared wavelengths. Thus the bilinear 
interpolation is suitable for the rural case and in this studv an acceptable method of 


interpolation for the marine case. 


TABLE? 
BILINEAR INTERPOLATION VALIDIT HE N 


ERROR(%)‏ )2 )ده( 
3.75am 121.008 3.75um 11. OOum‏ 
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Fig. 3.7 Monotonic behavior of the average mode | radius for marine and rural 
models from Table 1. 
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IV. RESULTS AND DISCUSSION 


After developing a technique to estimate extinction values in the infrared 
wavelengths from values in the visible and near-infrared, it is necessary to test its limits 
and sensitivity to several parameters that are included in the model. Five tests were 
designed with this purpose. The first test is a basic analysis of the technique itself in 
the retrieval of the particle size distribution. Two of the tests study how variations in 
the input parameters, extinction coefficient and relative humidity affect the estimated 
value. The fourth test looks at how deviations from the model size distributions affect 
the estimated extinction values. The last test compares how an extinction obtained 
from a three-mode particle size distribution compares to an extinction calculated from 


a model which assumes a two-mode particle size distribution. 


A. BASIC RETRIEVAL ANALYSIS 

The first test examines the error associated with the retrieval of a “perfect” 
satellite estimate of extinction after use of the interpolating technique. The procedure 
consists of the following steps: 


Í. ee an initial extinction values at 0.63 and 0.56 um. for a given TN, N2 
an 


ty 


Input extinction values into the first portion of the model to retrieve TN and 
A2 Which describe the particle size distribution. 


G3 


input TN and 2 into, the second section of the model to obtain an estimate of 
dun values for all four wavelengths of interest: 0.65, 0.36. 3.73 and 
.00 um 


4. Calculate the error from: 
error = 100(Pext estimated - Bext initial); ext initial (4.1) 


3. The procedure is done for both the marine and rural cases. 
The basic retrieval test for the marine case was run with values of TN and N2 set 


to 3,500 cm” 


and 0.035 respectively. In the rural case TN was set to 12,000 em’? and 
N2 to 0.000175. Humidity values of 25, 50, 75, 85, 96, 98 and 99% were used. 
Relative humidity values of 25, 75 and 85% needed interpolation by the model. As 
seen from the results in Table 3 the error associated with the retrieval technique is less 


than 2.5% for relative humidity values less than 99%. However, this value increases as 
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tne number of parameters that need to be interpolated increases. With all parameters 
interpolated, the error is still less than 10%. A larger error at the higher humidities 
was expected because this is where extinction changes the most as shown in Figs. 3.9 
and 3.10. The larger error at a relative humidity value of 25% is probably due to the 
fact that the interpolation for this value occurs over a large range, between 0 and 50% 
whereas at the higher humidities the difference is 10% or less. The error associated 
with a relative humidity value of 50% was the lowest and uses an uninterpolated 
matrix at a relative humidity where there is very little variation in extinction. This 
shows that the error due to interpolation could be lowered bv decreasing the range 
over which the interpolation is to be performed. This can be accomplished bv 


expanding the number of relative humidity values that are used by the model. 


TABLES 
ERROR ASSOCIATED WITH BASIC REINIEVALIEC MOM 


a. MARINE MODEL 


RH( %) RETRIEVAL ERROR(%) 
3.217 TOI 

25 5 2 3 

78 2 -0 3 

96 eM OM 

98 1S 1.9 

b. RURAL MODEL 

RH( % RETRIEVAL ERROR( % 
| 3. 75um 11. 00um 

| 0 0. 6 oe 

50 o OT 

85 S 1.5 

6.5 6.5 وو 


B.  EXTINCTION ERROR ANALYSIS 
The second test evaluated the sensitivity of the technique to errors in the 
extinction values used as inputs. Several assumptions were made, as discussed in 


Chapter II, for obtaining the initial input extinction value from satellite radiance 


measurements. Therefore, the extinction coefficient estimates may have significant 
errors. The test consisted of assuming the initiai extinction value is in error by £5, 10, 
20%. The procedure was as follows: 

For a given TN, N2, and RH, calculate Boyt at the four wavelengths of interest. 


2. Using the value obtained in step (1) calculate a value at 0.63 and 0.86 um 
S D Dy 主 + 50%. i ٢ Bext ۲ 


3. KRA above value into the first section of the model to obtain a retrieved TN 
an 


4. Input the retrieved TN, and N2 value into the second porna of the model to 
obtain an estimated extinction value at 3.75 and 11.00 utm 


Calculate the error as discussed in IV. A. 
Repeat procedure for 4 10 and 4 20%. 
The marine model was run for a TN value of 3,500 cm”? 
0.035. For the rural case TN was set at 15,500 cm? and N2 at 0.000135, using relative 


humidities of 25, 75, 96 and 98%. The extinction values were then varied Dv E10, 


and an N2 value of 


20 and 25%, as explained above for each humidity value. Table 4 shows the error 
associated with a variation in RH on an estimated extinction coefficient value at near- 
infrared and infrared wavelengths. As seen from the results, the magnitude of the error 
is basically determined from the accuracy of the extinction value within a couple of 
percentage points. This was expected since the retrieval method itself 1s accurate to 
better than 2% for most relative humidities. Slightly larger errors (3-59) for the rural 
model are observed for a humidity value of 98% where the value of extinction grows 
exponentially as seen in Fig. 3.10. The larger errors for the rural model at RH = 98% 
are observed for extinction errors of 20% or larger. This follows the inherent error of 
the retrieval technique at this humidity. Thus error in the predicted extinction is less 
than 3% for all humidities for both models after considering the initial extinction error 


and the error inherent in the retrieval teclinique. 


C. RELATIVE HUMIDITY ERROR ANALYSIS © 

Several assumptions were made in the development of the model for estimating 
extinction coefficients relating to relative humidity. The first assumption deals with the 
fact that relative humidity varies in the marine atmospheric boundary layer (MABL) 
with height. Relative humidity generally increases from the surface to the top of the 
MABL then decreases drastically to near zero at the top of the inversion layer as 
previously shown in Fig. 3.1. This model assumes input of an average value of relative 


humidity. The second assumption is that the relative humidity value is known. [t may 
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_____be_obtained from in situ measurements or in the future, perhaps, from satellite data 
(Kren, 1957). Changes in the relative humidity change the shape of the particle size 
distribution as shown in Fig. 2.3. In this study it was assumed that estimates of 
relative humidity value can be obtained within a certain degree of accuracy (2:599). A 
series of commercial humidity sensors were tested under laboratory conditions (Muller 
and Beekman, 1987). Test results show that for some of the sensors accuracies of 95% 
or better were achieved. Therefore, the test performed was in keeping with the 
accuracy reported for humidity sensors. 

The procedure for this study consisted of the following: 


l. Fora given TN, N2 and RH, calculate the initial extinction; for example enter 
an ‘observed’ (the real RH value) RH = 65%. 


tJ 


Use the extinction values obtained for 0.63 and 0.86 [im and input into the first 
section of the model. Enter an ‘assumed’ RH value that is H 5% (for example 
enter RH = 70%) to obtain TN and N2. 


N BPnier TN and IN2 irom above and the same RH value used in step 2 into the 
second portion of the model to obtain the estimated extinction. Following our 
DI ter an ‘assumed’ (due to errors in the instrument that measures RIT) 


/0. 


4. Calculate error as discussed in IV. A. 


Cn 


Se for -5%. Also repeat for different RH values for both marine and rural 
models. 


mie iaarinie case was run (or TN set to! 3200 cm * and N2 set to 0.0175 whereas 
the rural case values were 15.000 cm"? and 0.000135 respectively. Table 5 shows that. 
for the marine model, a 5% error in relative humidity values give rise to a range of 


5-149^ errors in Baq, at the near-infrared and infrared wavelengths. In the rural model 


ext 


a similar error in RH gives rise to estimated Baq, values being in error from 9-31%. As 


61 
seen in Table 5, the rural model is more sensitive to errors in relative humiditv at 11.00 
um wavelength than the marine model. Recall that extinction coefficients are a 
function of the particle size distribution. At 0.65 and 0.36 um the satellite data contain 
more information regarding the first mode in the rural case in contrast to more 
information on the second mode for the marine case. This can be seen if a vertical line 
is drawn for radii of 0.63 and 0.86 um on Figs. 4.1 and 4.2 until they intersect the solid 
line. The slope of the line at the point of intersection for the marine case is more 
dependent on changes in the shape of the second mode. Thus, a decrease in the slope 
infers a larger number of particles at the larger radii. As relative humidity increases the 
number of particles at the larger radii increases. At the higher wavelengths these larger 


particles contribute more to the extinction coefficient since they are more heavily 
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weighted. Since the information in Bag, about the second mode is limited in the rural 


ext 
model, a larger error is expected as observed in Table S. 


D. PARTICLE SIZE DISTRIBUTION ERROR ANALYSIS 
The next study investigates the error associated with errors in the particle size 
distribution obtained by changing the average radius of the second mode. The particle 
size distribution can be altered by changes in the relative humidity, as was shown in 
Ill. C. The distribution can also be changed by changing the average radius of the 
particles for the second mode, r2, of Eqn. 2.3. The second mode of the marine case 
was the only one tested since it is known that it consists of sea spray particles, which 
are more likely to varv in size due to variations in the wind. Changing r2 by. z: 2095, 
as seen on Fig. 4.3, moves the center of the second mode. The model was run with TN 
set to 4,000 cm’? and N2 set to 0.02. The procedure is as follows: 
l. Change r2 vaiue for a given relative humiditv value in the prediction model, bv 
35%. Calculate the initial extinction values representing a non-ideal particle 
size distribution. Reset the parameters. 


2. Enter extinction values into the first portion of the model to obtain values for 
TN and N2. 


C3 


Enter TN and N2 into the second section of the model to obtain the predicted 
extinction values. 


4. Calculate error as discussed in IV. A. 


Can 


Repeat for 410% and + 20% r2 variations in the marine case only. 
As we increase r2 we. change the distribution by increasing the number of 
particles with a larger radii. Table 6 shows that as we increase the variation in r2 the 
error associated with the positive deviations at }=3.75 jim increases from 0.5 to 10% 
while the error associated with the negative deviations increased from 6 to 25%. For 
入 = 11.00 itm the error associated with the positive deviations increased from 0 to 20% 
whereas for negative deviations the error increased from 9 to 41%. he error is 
greatest for negative deviations at the 11.00 utm wavelength. Again, at the 0.63 and 
0.86 um wavelengths the extinction efficiency function weighs the larger particles more 
heavilv for higher wavelengths as shown in Figs. 2.1 and 2.2. If a vertical line is drawn 
to the solid line for A = 0.63 and 0.86 jtm, as seen in Fig. 4.3, it can be shown that 


more information is contained in D,., at 0.63 and 0.86 um about the second mode for 


ext 
positive deviations. The net result is that the error is larger for the -20% r2 change. 
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TABLE 5 | 
RELATIVE HUMIDITY SENSITIVITY STUDY 


a. MARINE MODEL 


RH( %) PREDICTED ERROR(%) ۱ 
OBS ASSUMED 3. 75um 11. 00um 
65 70 11.8 9.0 
75 70 -10.3 -11.8 
85 90 8.6 12.6 
95 90 — -13.5 
b. RURAL MODEL 

RH( %) PREDICTED ERROR( %) 

OBS ASSUMED 3. 75pm 1 1 . 0 
55 70 21- 16.2 
75 70 1782 2E] 
85 90 Ez sid 
95 90 28.6 30.8 


E.  THREE-MODE PARTICLE SIZE DISTRIBUTION ERROR ANALYSIS 

The model for the estimation of extinction coefficients assumes that the marine 
particle distribution arises from particles of a continental and marine origin. In 
conditions of high winds this is altered by the introduction of larger size particles into 
the marine boundary layer. As the wind speed increases above 9 m/s, waves crest and 
particles are ejected into the atmosphere (Monahan er al., 1983). These particles have 
an average radius of 2 um and comprise the third mode of the particle size distribution 
(Gathman, 1983). The purpose of this test was to examine the effect of input of an 
extinction coefficient obtained from a variable wind regime and estimate extinction 
values using the assumed two-mode model. The model for the estimation of extinction 
coefficient was modified to allow addition of a third wind-dependent mode represented 
by variable A3 as shown previously in Fig. 2.4. The values of A3 equal to 1.0, 0.5 and 
0.1 were taken to represent the production of 10 jm particles at wind speeds of 
approximately 31, 27, and 15 m/s respectively. The shape of A3 is described bv: A3 = 
lOexp(.06 x wind speed - 2.8) and is a modification of the Navy Aerosol Model 
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TABLE 6 
PARTICLE SIZE DISTRIBUTION SENSITIVITY STUDY 


MARINE MODEL 


RH( 27) RZK 7) PREDICTED ERROR(%) 
SIE 11. OOMm 
30 TRO O5 oM O. O0 Lae 
SO SEO - 4.1 iem. AS 24.5 
S L-O -10.5 25.68 -20.6 40.5 
70 TAS QT 8.4 اس سح‎ Des 
70 EET ING - 4.0 c3 - 9.0 22 7 
70 E20 -10.0 Zr 219.0 oS 
95 S jos 6.4 "ic za 
85 310 - 1.1 د1‎ 5 - 6.1 19.4 
85 $20 - 8.2 107 "I 14.0 
90 =a نج‎ -- 8.8 -- ES 
90 10 -- 14.3 = 21 
20 T20 MAG 21.5 - 4.8 “6.5 


(Gathman, 1983). The model was run with TN set at 4,000 cm? and N2 set at 0.02 
for relative humidity values of 70 and 90% and TN was set at 4,135 cm? and N2 at 
0.0135 for relative humidity values of 60 and 8595. Different parameters were selected 
to test the sensitivity when matrices were not interpolated. The procedure was as 
follows: 


l. Obtain an initial extinction value of extinction for a given TN, N2, RH and A3 
value. For the first run set A3 to 1.0. 


2. Enter the extinction coefficients into the first section of the prediction model to 
retrieve a value for ITN and ۰ 


3. Enter TN, N2, RH and an A3 value of 0.0, which reverts the program back to a 
two-mode distribution, into the second portion of the prediction model to 
obtain the predicted extinction values at 3.75 and 11.00 pm. 


Calculate the error as discussed in IV. A. 


Repeat for A3 values of 0.5 and 0.1. Also repeat for several relative humidity 
values for the marine model only. 


‘From H. Hughes by, personal communication, present affiliation: Ocean and 
Be cap Sciences Division, Naval Ocean Systems Center, San Diego, CA 


e£ 


47 


Wirin the addition of a third mode we observe in Table 7 tnat the error 's ¢reatest 
at the nigher wind speed, up to 82°% for A3=1 and X=11.00 pm and up to 39% at 
A= 3.75 pm. Results show that the errors decreased as relative humidity values 
increase. ‘the addition of more particles at the higher radii for lower RH values 
influences the result more and thus the higher error at the lower RH values. Also, 
since the higher radii are weighted more by the extinction efficiency function for the 
higher wavelength a higher error is observed at 11.00 um. There is a slight increase in 
error at 60% RH due to the added interpolation: step over a 20% range of relative 
humidity. The interpolation error is not observed at 35% RH since the interpolation is 
over a 10% spread and the error associated with the introduction of the third mode 


dominates. 


TABLE 7 
EFFECT OF A THREE-MODE PARTICLE sl4E ۱ لرل‎ 


MARINE MODEL 


RH( % ) A3 PREDICTED ERROR) 
3, ٢ 1 1. ۰ 
60 1:6 =S 950 -81.7 
60 3 ۰30. 2 -75.8 
60 0 1 - 7.1 -47.2 
70 0 -30.2 EE 
70 0,5 -—22 9 -67.0 
70 0 - 4.0 sa 
85 TO -- -- 
85 8 د‎ = Ls =52 د‎ 
So Ol 8. ۵ =: 3.5 
90 || = ال د‎ -34.5 
90 0 = Ses -22.4 
90 01 so ES 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A method for estimating extinction coefficients in' the near-infrared and infrared 
wavelengths using channels 1 and 2 of the NOAA-7 AVHRR was developed. The 
model consists of two sections. Input of the type of aerosol model (marine or rural) 
and the effective relative humidity value for the boundary layer defines the average 
mode radius. Input of the extinction coefficient at 0.63 and 0.86 um determine the 
magnitude of the retrieved parameters, TN and N2, necessary to describe the 
distribution of aerosols. The second section uses these values and the effective relative 
humidity value for the boundary layer to estimate the extinction coefficients at 3.75 
and 11.00 um. 

Five tests were developed to determine the limits and sensitivitv of the model. 
The first test examined the error associated with the retrieval of a ‘perfect’ satellite 
estimate of extinction. Given an initial value of extinction at 0.63 and 0.86 um for a 
given TN and N2, results show that the retrieval error is the smallest at the middle 
humidity values where the extinction coefficient, as a function of relative humidity, 
changes the least and interpolation occurs over a relatively small range (10% spread). 

Results for the second test show that the models’ accuracy is basically determined 
from the accuracy of the input extinction coefficient values. Error in the predicted 
extinction is less than 3% after considering the error in input extinction for all 
hunudities tested and for both models. 

Relative humidity affects the distribution of large radii particles in the marine 
atmospheric boundary layer. A 5% error in the value of the relative humiditv for the 
area of interest will lead to errors as high as 14% for the marine model and 31% for 
the rural model in the estimated extinction coefficient values at the longer wavelengths. 
This test shows that errors in relative humidity affect the rural model the most, since 
the information contained in the second mode is more limited for the rural model for a 
given relative humidity. 

The particle size distribution can also be changed by increasing the average size 
of the particles in the boundary layer. Thus changes in the mode radius by 5, 10, 20% 
will result in an error in the estimated extinction coefficient value of up to 12, 24 and 
40% respectively. The highest errors are for calculations at the 11.00 um wavelengths 


since the extinction efficiency function weighs the larger radii the most. 
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«he last test evaluates ine performance of a two-mode particie size distribution 
mode! azainst the Navy aerosol three-mode model. With the addition of a third mode, 
we observe that the error decreases as the wind speed decreases. Since the third mode 
contains very large particles, the error is highest for estimates at 11.00. jtm as expected. 

Three general conclusions can be made from this study. First, estimates of 
extinction at infrared wavelengths rely heavilv on the accuracy of the input extinction 
coefficient values at 0.63 and 0.86 um. Therefore, care must be exercised in obtaining 
the best possible value of extinction from satellite estimates of optical depth and 
NOLAPS estimate of boundary layer depth. The second is that extinction 1s dependent 
on particle size distribution and thus the shape of the distribution used in the model 
will affect the estimated value directly. And finally, factors affecting the particle size 
distribution, such as relative hunudity and wind speed, must be taken into account. 

It must be noted here that the model as it stands now is limited to certain TN 
and N2 values. To be used operationally these parameters and the associated 
extinction matrices would have to be expanded. Error in the interpolating technique 
can be decreased by increasing the number of relative humidity values that are used in 
the model. Also, due to large particles being generated by high winds. the model 
should be modified to allow for variable winds similar to the Navv aerosol model. 
Lastly, methods for estimation of relative humidity, to accuracies of 95% or better, 
need to be developed for meteorological satellites to facilitate measurement of relative 


humidity in remote areas and thus provide better input into the model. 
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